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Abstract 

We propose a solution to the fJ-fBfj. problem in gauge-mediated SUSY breaking, 
which does not suffer from the SUSY CP problem and is consistent with the solution 
to the strong CP problem. The model is based on Zs-invariant NMSSM with addi- 
tional vector-like matter charged under a hidden gauge group as well as the standard 
model gauge groups. The dynamical scale of the hidden gauge symmetry is set to 
be around 10 GeV. We show that this simple extension of the NMSSM resolves the 
domain wall problem and the f^/Bfi problem without a dangerous CP angle. The 
relative sign among gaugino masses and the n parameter can be preferable in terms 
of the muon anomalous magnetic moment. We also discuss cosmological issues of 
this model, especially the effects of long-lived particles in the hidden gauge sector. 
The hidden glueball may cause a late-time entropy production, which opens up a 
possibility that the light gravitino can be a dominant component of dark matter 
while leptogenesis scenarios work successfully. There is a region where the elec- 
troweak symmetry is successfully broken with ~ 10 TeV stops, and the Higgs mass 
of ~ 125 GeV can be explained. 



1 Introduction 



Gauge-mediated SUSY breaking (GMSB) is the one of the attractive mechanisms to 
transmit the SUSY breaking effects from the hidden sector to the visible sector jl]. In 
GMSB models there are no sources of CP violation and flavor violation in sfermion and 
gaugino sectors, since SUSY breaking effects are purely transmitted by gauge interactions. 
However, there are difficulties for generating viable /i term and Bfi term. The first is the 
so-called ^/B^ problem, which states that it is difficult to generate correct size of both /i 
and B parameters for the electroweak symmetry breaking in GMSB models. The second 
difficulty is that there generally exists a CP violating phase among /i term, Bfi term and 
gaugino masses. Unless the relative CP phase is smaller than O{10~^), it induces too 
large electric dipole moments of the electron and hadrons, which are excluded by the 
experiments [2]. 

There is an elegant solution to the fi/B^ problem in the framework of Next-to-Minimal 
Supersymmetric Standard Model (NMSSM) [3]. By introducing additional vector-like 
matter charged under SM gauge groups to NMSSM, viable and Bfi are generated with- 
out introducing an additional CP violating phase [H El [6l [TJ [8]. It also eliminates the 
cosmological domain wall problem, which generally exists due to the spontaneous break- 
down of the symmetry in NMSSM, because the symmetry is explicitly broken by the 
quantum anomaly induced by the vector-like matter p] . However, this solution is not com- 
patible with the Peccei-Quinn (PQ) symmetry for solving the strong CP-problem pUf [TT]. 
If this scenario is combined with the PQ mechanism, a linear combination of the original 
Z3 symmetry and the PQ symmetry leads to another non- anomalous Z3 symmetry, and 
the domain wall problem is restored [9l |8]. 

In Ref. [H] we constructed a model based on NMSSM with additional vector-like matter 
charged under a hidden gauge group and showed that such a model induces sizable fi/Bfi- 
term, solves domain wall problem and is consistent with the PQ mechanism. This model, 
however, was rather involved : it contained a hidden messenger sector, which mediated 
the SUSY breaking effect into the hidden matter sector. 

In this paper we propose a simpler extension of the NMSSM in GMSB. We only 
introduce additional vector-like matter charged under both the SM and hidden gauge 



2 



groups. The messenger sector is as minimal as the the original GMSB model. In this 
setup, we will show the fJ^/Bfi problem, the SUSY CP problem and the domain wall 
problem are all solved, consistently with the PQ mechanism for solving the strong CP 
problem. Moreover, the relative phase of the /x-term and the gaugino masses can take a 
favored sign from the viewpoint of the SUSY explanation of the deviation of the muon 
g — 2 (anomalous magnetic moment). 

The novel difference from the previous model [8] is that cosmology of the present 
model may be significantly modified from the standard one because of the existence of 
the long-lived particles charged under both the SM and hidden gauge groups. Moreover, 
the hidden glueball may become a source of a late-time entropy production, which might 
be welcome in GMSB where the upper bound on the reheating temperature after inflation 
is stringent. The gravitino is a good dark matter candidate after the dilution by the decay 
of hidden glueball. 

This paper is organized as follows. In Sec. [2] our model and its basic structure is 
described. In Sec. |3] we describe phenomenology of the model, including the mass spec- 
trum, collider signatures, various constraints and implications on muon g — 2. In Sec. H] 
cosmological aspects of this model is discussed in detail. Sec. |5]is devoted to conclusions. 



In this section we introduce a model of NMSSM in GMSB with extra matter, which 
avoids the domain wall problem and is compatible with the PQ solution to the strong 
CP-problem. 



2 NMSSM in GMSB with extra matter 



2.1 Model 



We consider the following Z3 invariant superpotential: 



W = WnMSSM + W^GMSB + W^extra 



(1) 



where 



W^NMSSM 



XSHuHa + —S + WMSSM-Yukawa , 



(2) 
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Table 1: Charge assignments on chiral superfields fields in the model under the Z3, 
R-parity (+ : even, — : odd) and S\J{N)h- 5^ and 10m are the MSSM matter fields. 

W^GMSB = XM + W^hidW, (3) 
W^extra = ko' S D' D' + S L' L' . (4) 

The charge assignments of chiral superfields are shown in Table [1] VFnmssm is the super- 
potential of the NMSSM, where singlet Higgs, up-type Higgs and down-type Higgs are 
denoted by S, and H^, respectively. The effective yU-term is induced by the vacuum 
expectation value (VEV) of the singlet field, fies = A {S). Wgmsb is responsible for the 
gauge mediation of the SUSY breaking, where messenger superfields are denoted by 
and ^, which are 5 and 5* representation of SU(5) grand unified theory gauge group. 
The SUSY breaking field X is assumed to develop a VEV as X = Mmess + ^'^-Fx due to 
the dynamics induced by lyhid(X) (and Kahler potential). 

Finally, D' and L' {D' and L') are extra vector-like matter, which are charged under 
the SM gauge group as well as an S\J{Nh) hidden gauge group. They naturally induce 
the VEV of the 5* field through the effects of renormalization group evolution on the mass 
of S, and also play an important role to break the symmetry via anomaly. As will 
be discussed, the hidden gauge group S\J{Nh) is necessary to break the symmetry 
in the presence of the PQ mechanism. Note that Nh should be either 2 or 3, since the 
perturbative gauge coupling unification fails for larger Nh- As we will see, it turns out 
that the dynamical scale of the hidden gauge group. Ah, should be about 10 GeV. The 
mass spectrum of the model is discussed in Sec. |3l 

2.2 Domain wall and strong CP problem 

It is known that the Z^ invariant NMSSM suffers from the domain wall problem since the 
Z3 symmetry is spontaneously broken after the electroweak phase transition. One may 
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introduce a small explicit Z3 breaking term in the superpotential, but this in general leads 
to a tadpole term which tends to destabilze the gauge hierarchy [121 [13] Fl 

In the present model, however, the Z3 is anomalous at the quantum level due to the 
additional vector-like matter. By integrating out D' , D', L' and L', a Z3 breaking term 
is induced as 



where G'^^^ is the hidden gluon field strength and 6h is the sum of the strong phase 
of the hidden gauge group '&\]{Nh) and phases induced by the mass of extra-matter, 
which is in general expected to be 0(1), and counts the number of matter fields with 
fundamental and anti-fundamental representation of SU(A'^i^) (A^ = 5 in our model). We 
denote a singlet CP-odd Higgs by a^, and vs = (S) is the VEV of 5*. Apparently, this 
term violates the Z3 symmetry while the rest of the Lagrangian is invariant. After the 
hidden QCD phase transition, there arises a potential for as which explicitly breaks 
symmetry [T7] 



where f{x) is a periodic function which satisfies f{x + 27r) = f{x) and Ah is the strong 
scale of the hidden gauge interaction [18]. Therefore, by comparing it with the tree level 
Za-invariant scalar potential, we see that the degeneracy of the three distinct vacua are 
broken if an integer A^ and 3 are relatively prime. The bias among the potential energies 
of the original three vacua are AV ~ A%. As exphcitly shown in Ref. [8], A^ > C(l) MeV 
is sufficient to make domain walls unstable so that they decay well before dominating the 
Universe. Note that the potential (|6]) slightly shifts the position of the minimum of as 
from zero, which induces a small CP violation. We will discuss it in Sec. 13.51 

Here let us see that our model is compatible with the PQ solution to the strong CP 
problem [TOllIT]. Since the MSSM sector does not exhibit a PQ symmetry, we introduce 
a new sector in the superpotential. 




(5) 



V^,r^A%f{9H + Nas/vs), 



(6) 



W = k^pqQpQQpQ, 



(7) 



^ See e.g., [131 [TH [Kl [IS] for models without the tadpole problem. 
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where $pq is a gauge singlet and Qpq (Qpq) are (anti-)fundamental representations of 
SU(5) 



The global U(1)pq charges are assigned as $pq(+2), Qpq(— 1), Qpq( — 1). After 
<I>PQ develops a VEV of fa ~ C(10^° - 10^^) GeV, the U(1)pq is spontaneously broken 
and then an almost massless mode corresponding to the angular component of $pq, axion 
(apq), appears. In this paper we do not specify the mechanism of the stabilization of the 
PQ scalar along the flat direction associated with the U(1)pq symmetry. Although it is 
irrelevant for the following discussions, we here show one example. A simple way is to 
introduce a superpotential of the form W = $pq$pq/M"~^ with a cutoff scale M where 
l>PQ has a PQ charge -2n. Then the PQ scalar is stabilized at ($pq) ~ (mpQM"''^)^/^""^) 
and (5>pq) = 0, where mpq is the SUSY breaking mass for the PQ scalar. By choosing n 
appropriately, we can obtain a desired PQ scale = ($pq). 

The U(l)pQ is anomalous under QCD. Since the Z3 is also anomalous under QCD, 
the instanton effects induce the following potential, 

VpQ ~ A^cD [1 - cos (apq//, + Nnas/vs + 9o)] , (8) 

where 60 is the bare strong CP angle. Independently of the value of as, the axion always 
dynamically cancels the CP angle and resolves the strong CP problem. Notice that, 
without a Zs-breaking potential induced by hidden gauge group (Q, we could always 
choose apQ so that Vpq = for each Z3 vacuum, i.e., there would remain an unbroken 
Z3 symmetry. This is why we need a strong hidden gauge symmetry to make the 
anomalous. 

3 Mass spectrum and phenomenology 

Let us now discuss the mass spectrum and phenomenology of the modeljfl The model is 
parameterized by the following parameters: 

Fx 

Mmess, A = — , A, K, ko', ki', Nh, Qh- (9) 



^ We introduce only one set of such matter pairs. If more matter is introduced, the gauge couphng may 
blow up below the GUT scale. Note also that color anomaly number is equal to one for the model with 
only one pair, and hence there is no axionic domain wall problem even if the PQ symmetry is restored 
after inflation [19j . 

^ The analyses in this section are out of date in light of the recent discovery of the 125GeV Higgs 
boson after the submission of this paper. See "Note Added" for updated analyses. 
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In our analysis, we fix k and koi in terms of tan/3 and the Higgs VEV, similar to the 
case of MSSM where /i and Bfi are fixed. As described below, the mass squared of the 
gauge singlet, m|, which should be negative for successful EWSB, is controlled by the 
parameter kr)/. The parameter ki' does not affect the EWSB as significantly as kc,/ does. 
Therefore in numerical calculation, we fix it to be a small value, since otherwise there 
appears an unbounded vacuum, as described in Sec. 13. 3[ should be 2 or 3 in order to 
keep the perturbative coupling unification. As will be discussed, viable phenomenology 
and cosmology require Ah ~ 10 GeV, which corresponds to (^//(ITeV) ~ 1.2 (1.4) for 
Nh = 3 (2). Thus, there are essentially 4 parameters left, Mmess, A, A, and tan/3. 
Below the messenger scale, the following soft terms are induced, 

-£soft = ml\S\'^ + mjjjHul'^ + mjjjHdl'^ 

+ml,\Df + ml,\Df + ml,\Lf + ml\L'\^ 
+{AxXSHuHd + A^KSy3 + h.c.) 

+ {AD'kD'SD'D' + AL,kL,SL'L' + h.c). (10) 

where the soft terms of the squarks, sleptons, and gauginos are omitted for simplicity. 
The soft masses for the extra matter, m|,,, m^,, roj^, are generated in a similar way 
to the MSSM squarks and sleptons via GMSB mechanism, since they are charged under 
the SM gauge groups. In the presence of this extra matter, the beta-function for m| 
receives additional contributions, 

(Stt^)^^ 3 3NHkl,{ml,+m%+ml + \AD'\^) 

+2NHkl,{ml, + ml, + m| + \Al'\'^) (11) 

where Q is the renormalization scale. This induces a negative m| and the VEV of S field 
accordingly, thereby solving the ^/B^ problem. Since m|,, and m^, are larger than m|, 
and m|, in GMSB, the m| and (S) are mainly controlled by the parameter ko'- 



In the numerical calculations, we have used the IN 
to include the effects of additional vector-like fields 



MSSM Tools [20J, which is modified 
and to include two-loop effects of 



* ku' is determined iteratively so that the conditions for the EWSB are satisfied with predicted soft 
mass parameters. 
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PI I 267.9 I 4559.9 I 571.4 I 617.6 I 4411.0 I 5107.2 



P2 332.4 4966.9 561.8 619.9 4794.2 5515.0 
P3 402.4 " 4471.9 luT 771.8 4460.1 "5198^ 

Table 2: The mass spectra of some model points are shown. The input parameters are A, 
^mess, A, tan/3, and Nh. The parameters, k^)', k and /iefr are determined by iteration. 
kii is taken as k^i = 0.002 for all numerical calculations. The fermion masses of L' and 
D' are denoted by m^^, and m^^, respectively. All masses are written in units of GeV. 



additional hidden gauge group (see Appendix A). The mass spectrum, the LEP constraint 
and the branching ratio of the inclusive B meson decay are calculated in NMSSM Tools. 
Now let us see various constraints and implications of the model. 

3.1 NMSSM mass spectrum 

Typical viable parameter regions are shown in Fig. 1-4 and the mass spectrum of some 
model points are shown in Table El Various constraints and phenomenological implications 
of the model will be discussed in more detail in the following subsections. In the viable 
region we are interested in, the SUSY breaking scale is Fx ^ 10^^ GeV^ and hence the 
gravitino is heavier than 0{1) keV. The lightest SUSY particle within the NMSSM sector 
is the lightest neutralino, which is stable inside the detector. In these cases, SUSY searches 
limit the squark masses, as rriq > 1100 GeV for the gluino mass rrig ~ 900GeV, and 
rriq > 1200 GeV for rrig ~ 800GeV [2T1I22]- Figs. Eland H] show contours of the gluino mass 
and squark mass, respectively. There are viable regions where the constraints from the 
SUSY searches are avoided while the deviation of the muon g — 2 from the SM prediction 
is within la level (see also PI, P2 and P3 in Table Note that for Nh = 3, the 



8 



one-loop renormalization group coefficient of the gluino mass vanishes and the two-loop 
contribution is important for deriving the mass spectrum. (See Appendix A). 

In typical viable points, our model predicts that the lightest CP-even Higgs, hi is 
singlet like and is as light as (9(10) GeV|f| Therefore it is excluded by the LEP results if 
its mixing with the SM-like Higgs is too large. Since larger A leads to larger mixings with 
Hu and if^, there is an upper bound on A to avoid the LEP constraints on the mass of hi 
and the coupling, Z — Z — hi [23]. In Fig. |2l we show the constraint on A on A-A plane. In 
order to avoid the LEP constraint, A should be as small as (9(10^^). On the other hand, 
the second lightest CP-even Hi^gs is SM-like, and its mass about 117GeV by including 
top/stop radiative correctionsjj 

In most of the viable regions, the lightest CP-odd Higgs is singlet like and is as light 
as a few GeV. The lightness originates from the smallness of k (~ 10~^) and A^, which 
corresponds to an approximate PQ symmetry^ Smallness of k is guaranteed from the 
requirement that the correct electroweak symmetry breaking is obtained when A ^ 1 (i.e. 
{S) ^ ITeV) and the soft mass parameters, |?^//^| and |^«||, are not much larger than 
10^ GeV^. Such a very light CP-odd Higgs can affect the prediction of Br(_B — )• X^/i^/i^) 
significantly, as described Sec. 13.41 

3.2 Anomalous Magnetic Moment of the Muon 

The experimental value of the muon g — 2 deviates from the SM prediction with about 
3a level [25l [26] . In Ref . [25] it is calculated as 



(a^)EXP - (a^)sM = (26.1 ± 8.0) ■ 10"^° (12) 

where the subscripts EXP and SM refer to the experimental value and the SM prediction. 
Interestingly, the deviation can be naturally explained by SUSY contributions when tan /3 



^ The solution with /icff < predicts heavier Higgs masses (see Refs. [IllH]). 

^ We have checked that FeynHiggs [53] also predicts the SM-like Higgs boson mass of ~ 116.5 GeV, 
with the parameters in PI and P2. 

^ This is different from the PQ symmetry discussed in Sec. 12.21 for the solution to the strong CP 
problem. 
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Figure 1: The contours of g-2 and allowed regions of the parameters. A=100 TeV, 120TeV 
and 140TeV for (a)(b), (c)(d) and (e)(f), respectively. In panels (a)(c)(e), A = 0.006 and 
Nh = 2, and in panels (b)(d)(f), A = 0.005 and Nh = 3. In orange (yellow) regions, 
muon g-2 is explained at la {2a) level. Pink regions are excluded by the LEP constraint. 
Blue regions are excluded by i? — )• Xgfi^fi^. Green regions are excluded by unbounded 
vacuum. Gray regions are not consistent with the experimental data of the muon g-2, 
since there is no solution with /lefj > 0. The contours are drawn with gH{fnq) = 1.2, where 
Trig is a squark mass. ^ „ 
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Figure 2: The allowed region of A is shown. In the left panel, we take = 2, tan /3 = 40 
and Mjness = lO^GeV and in the right panel, Nh = 3, tan /3 = 40 and Mmess = lO'^GeV. 
gni^q) is taken as guijnq) = 1.2. The meanings of the colored regions are same as those 
in Fig. 1. 




Figure 3: The contours of gluino mass are shown. The numbers are expressed in the 
unit of GeV. In the left panel, A'^^^ = 2, A = 0.006 and tan /3 = 40 and in the right panel, 
Nh = 3, a = 0.005 and tan/3 = 40. gn is taken as gnirrig) = 1.2. The meanings of the 
colored regions are same as those in Fig. 1. 



11 



10 



10 



10- 



> 
o 



J 10« 





1200 












\ 

1400 


\1600 








' r 


























1 










10' 

80 90 100 110 120 130 140 80 90 100 110 120 130 140 150 

A (TeV) A (TeV) 

Figure 4: The contours of squark mass are shown. The numbers are expressed in the 
unit of GeV. The parameters and colored regions are same as in Fig. [3] 



is (9(10) [27J. The SUSY contribution to the muon — 2 is approximately given by [28j^ 



9l 



m. 



32vr2 ml,, 



tan/3- 



M2/i| 



(13) 



It is enhanced for large tan/3 and suppressed for large ml,,. In order to explain the muon 
g-2 anomaly, the sparticles (sleptons) should be sufficiently light. Note that the relative 
sign between M2 and fi is also important. For M2/i < 0, SUSY contributions tend to make 

more discrepant from the SM prediction. In Fig. 1-4, contours of g — 2 are shown. In 
the orange (yellow) regions, the muon g — 2 agrees with the experimental result within 
Icr [2a) level. The viable region consistent with the muon g — 2 corresponds to Fx — 
IO12 _ 10^6 GeVl Therefore, the gravitino mass is mg/a ^ Fx/Mp ~ C(0.001 - 10) MeV. 
Coincidentally, such a light gravitino is also preferable from the viewpoint of cosmology 
as described in Sec. HI 

Since SUSY searches excluded the light sparticle mass regions in GMSB, a large tan/3 
is necessary in order to explain the muon g — 2. For large tan/3, only viable region with 



^ This approximation is valid only when the diagram in which and x dominantly contribute to 
the loop among other SUSY contributions. 
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successful electroweak symmetry breaking is |k/A| ^ 1 and A ^ 1 as long as soft masses 
are induced by GMSB effects [7]. This is the so-called decoupling limit of the NMSSM. 
However, very small values of A are also constrained from the vacuum structure, as we 
will see below. 

3.3 Vacuum structure 

In the NMSSM, the vacuum structure is complex compared to that of MSSM, due to the 
presence of another singlet Higgs; there exist unwanted local (global) minima in general 
(c.f. [29]). Therefore the depth of these minima should be compared to that of the wanted 
minimum so that the successful electroweak symmetry breaking is realized. 

The relative depth of the wanted minimum, which we call "true" minimum hereafter, 
measured from the origin S = Hu = = is roughly given by, 

Vtrne ^ -(m|)V(4/s:2) ~ -10l^GeV^ (14) 

where we neglect small A-teim contributions. Since we consider the decoupling limit 
where vs is as large as O(IOO) TeV, the VEVs of Hu and are approximately neglected 
in the following analysis. Now let us see other local or global minima existing in the 
present model. 

3.3.1 Avoiding the color breaking minimum 

First note that there may be additional constraints on the vacuum structure once the extra 
matter is included |H]. The mos stringent constraint comes from requiring the absence 
of unbounded direction. The unbounded direction exists when the soft masses of D' and 
D' and/or kfj' are small. Along the direction of {S) = , | | = | {H^) \ = vh and 
I {D') I = I <^Z)') I = t>£i, the scalar potential reads 

V = (-lAI-y^ + \kD'\vl,^ + (m^^ + m]jjv]j + (m^, + m^,)^!)- (15) 
The first term vanishes for v'j-, = IX/k^ilvfj. Then the scalar potential becomes 

V = {m\^ + + (m^, + mr,:)\\lkD>\)v\. (16) 
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Apparently, large fc^/ leads to an unbounded vacuum, since the contributions from soft 
masses of D' and D' are suppressed, and negative m|^^ dominates over other contributions. 
Thus fc^)/ cannot be very large for the stability of the vacuum. Similarly, large k^i is also 
not allowed. Since m\i and m|, are smaller than m^, and the upper limit on fc^/ 
is more stringent. Therefore the SUSY mass of the extra lepton-like states, L' and L', 
should be as small as (9(100GeV). |^ In numerical calculation, we fix k^i = 0.002. 

3.3.2 Stability of the true vacuum 

There is a potentially dangerous local extremum along (H^) 7^ and (S) = (H^) = 0. 
The relative depth of the potential along this direction is estimated as 

Vh. = ~ -lO^'GeV. (17) 

92 

for —rn^j^ ~ |yUp ~ lO^GeV^. Thus the desired (electroweak symmetry breaking) mini- 
mum is deeper : Vewsb < Vh^ in the parameter space which we are interested in. 

Notice that there is a deeper minimum along the D-flat direction for the Higgs: 
I {H^) I = I {H^) I = vu and {S) ~ 0. Along this direction, we have the deepest min- 
imum, 

Vh ^ -{ml^ + m|J V(4A2) ~ -lO^^GeV^ (18) 

This is deeper than the desired minimum. Thus the desired minimum that we want is 
metastable and we need to take care about stability of the desired vacuum. 

First, let us see the stability of the desired vacuum. The desired minimum at vh^ <^ 
vh^ vs is separated from the global minimum by the potential barrier, which is of the 
order of Harrier ~ XAxv^ ~ lO^^GeV, where v ~ 100 TeV. Thus the Euclidean action for 
the bounce solution connecting the two minima is estimated to be 5*4 ~ lOVbarrier/'^l/^ ~ 
10^, which is sufficient for avoiding the decay of the metastable vacuum into the global 
minimum on a cosmological timescale [5U] . 

The next issue to be discussed is whether or not the Universe prefers the desired 

metastable minimum rather than the deeper global minimum. At high temperature regime 

^The breaking of the GUT relation between fc^/ and k]jf may come from the existence of a higher 
dimensional operator, e.g., S" 5(Z)^ L')S(24)5/Mp (S(24)_is the SU{b) GUT breaking Higgs); additional 
corrections, e.g., ~ 2 x 10~'^SD'D' and ~ — 3 x may arise. 
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of T ^ 1 TeV, all particles with masses less than ~TeV are thermalized and hence they 
give thermal masses for S, and Hd- This stabilizes these fields at the origin : S = = 
Hd = 0. As the temperature decreases, tachyonic directions appear. Thus the scalars roll 
down toward the direction which becomes tachyonic first. By comparing the thermal 
mass for the 5* direction, ms(T)^ ^ SNnkj^/T"^ and that for the HuHd direction, 
rriH^HaiTy ~ [3{yf + yl) + [l^g"^ + ^g'"^) /2]T'^ /12, we obtain the following condition for 
the 5* direction to become tachyonic first : 



I mil \mij^ + mH^ 



Nnkl, (^1 + ^2) + 5^2/2 + 5^/2/6' ^ ) 

where we have assumed D' , top and bottom quarks, gauge bosons, gauginos as well as 
higgsinos are in thermal bath. Actually this is satisfied for typical parameters in the 
model. Therefore we conclude that the Universe relaxes at the desired minimum toward 
vs ~ Tnsl ~ C(IOO) TeV with vh^ ~ ~ at the phase transition. The /i and 
terms are generated by the VEV of S, and the model effectively looks like the MSSM. 
After that, Higgs fields found the desired vacuum where ■\Jvjj~\-vj^^ = 174 GeV. This 
vacuum is stable against decay into the global minimum, as shown abovec^ 

3.4 Bs ji^ji^ and B Xsii^fi^ 

When the lightest CP-odd Higgs, Oi, is very light, there are sizable contributions to 
the processes involving B meson [311 ESI |33l |3l]. The flavor violating coupling among 
bottom, strange quarks and Oi is generated by non-holomorphic Yukawa couplings, which 
are induced by picking up the CKM matrix element and soft masses at one- loop level • 
Since both A and k are small in our model, the mixing between the MSSM-like CP- 
odd Higgs and the singlet like CP-odd Higgs is as small as (9(10~^ — 10~^). Therefore 
the contributions to — )■ fi~^fi~ mediated by ai are small, resulting in the prediction of 
Bt{Bs —7- fJ^^fi') ~ O{10~^). The branching ratio is about one order of magnitude smaller 
than the limit from the LHCb result, Br(S, ^ fi+fi-) < 1.5 x 10^^(95% C.L.) [36j. Ill 



Precisely speaking, the first order phase transition into the global minimum may slightly precede 
the second order one. This does not modify the above result much as long as the phase transition into 
the true minimum occurs much earlier than the second order phase transition into the global minimum; 
i.e., as long as the condition (fT9|) is well satisfied. 

Although, the upper bound for Br{Bs fx^ iJi~) is updated as Br(_Bs < 4.5 x 

10~^(95% C.L.) the most of the parameter space can satisfy this new bound. 
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The constraints from AM^ (g = s, d) are also loose with such a small mixing [3^ . 

On the other hand, the inclusive B decay, B — )■ X^fi^fi^ gives a more stringent con- 
straint and excludes some regions of parameter space, even though the mixing is very 
small. When the mass of the lightest CP-odd Higgs is in the ranges of 1 GeV < M^+^- < 
2.4 GeV and 3.8 GeV < M^+^- < 5 GeV, the ai, which mediates this process, can be on- 
shell. Therefore the regions with these CP-odd Higgs masses are excluded in most cases. 
The blue regions shown in Figs. 1-4 are excluded by this constraint. Thus the regions 
with small messenger scale and small tan /3 are not favored. Note that the region between 

2.4 GeV < M^+^- < 3.8 GeV can not be constrained, since it is difficult to estimate the 
branching ratio due to the charm quark resonances. (The corresponding regions are those 
between two blue strips in Fig. 1-4.) 

3.5 Small CP violation 

In NMSSM without CP violating parameters, the CP violating extrema in the potential 
is not minimal but maximal, and hence CP is conserved. It should be noted, however, 
that the potential ([6]) slightly shifts the position of the minimum of from zero. Since 
there is no SUSY CP phase only for as = 0, the inclusion of this term induces a small 
spontaneous CP violation. The induced CP angle is estimated as 

- ^ - 4^"- io-.„ (-A^y ('-^r (™)\ (20) 

{S) miys VlOGeVy \ TTia^ J \ vs J 

The phase of the singlet VEV, 56s, gives a phase to that of the effective yU term, /^efr = Af 5, 
which is constrained by EDM experiment [2]. Since the upper limit on the CP phase is 
0(10-3) - O(10~^) for sparticles of 1 TeV and tan/3 = 0(10), the predicted 59s is well 
below the experimental bound unless Kh is much larger than 0(10) GeV. 

3.6 Hidden gluon and glueball 

The hidden gauge sector contains no light matter. Thus it is expected that the hidden 
glueball is the lightest particle in the hidden gauge sector, and hence it can be long-lived 
on a collider (and even cosmological) time scale. Let us see properties of the hidden 
glueball. 
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Once a hidden glueball is formed, it can decay into SM particles through loops of extra 
matter |3H1 139] . For example, there are effective dimension 8 operators after integrating 
out the extra matter, 



ana' 

-L-eff ~ — 4 — 



F.^F^^'Ti {GupaG'^) + F;..F'^"Tr (Ghp.&h)] , (21) 



where Gh denotes the field strength of the hidden gluon and an is the coupling constant 
of the hidden gauge symmetry. The first term induces the scalar glueball decay and 
the second term the pseudo-scalar glueball decay. Both have a similar lifetime, which is 
estimated as 

.(,„^,„.6xl0-=.c(l^)°(i!^/, (22) 

where we have simply approximated all the mass scales appearing in the pure hidden 
gluon sector, the glueball mass and its decay constant, to be Ah- Thus Eq. ( l22l) should 
be regarded as only an order-of-magnitude estimation, but it is sufficient to see that the 
lifetime can be short enough to be free from the BBN constraint, and long enough to pass 
through the detector at collider experiments without leaving signals except for missing 
energies. Similarly, they can decay into SM gluons via dimension 8 operators, but in this 
case it is D' that mediates the decay, and it is heavier than L'. Due to the high powers 
of rriDi in the expression of lifetime, it may not be the dominant decay mode. 

Moreover, the one-loop coupling of hidden gluons to the singlet scalar leads to the 
following effective dimension 7 operators 



'^'^ //Tr {Gh,.G'h) + ^^/75/Tr {Gh,.G''h 

mi — Arr V 



(23) 

iTivs -Ah ' " - Aj^ \ ' " / j 

where / denotes SM fermions which are mediated by exchange of scalar and pseudo-scalar 
Higgs bosons. Here e^h and eaga^ represent mixings between the singlet (pseudo) scalar 
bosons and SM (pseudo) scalar Higgs bosons. For example, the lifetime of the scalar 
glueball into the charm pair is estimated as 

^ e.) 5 X 10-=ec (— ) (^) (^) (- ) , (24) 

and similar expression holds for the pseudo scalar glueball. This is close to that induced by 
dimension 8 operators (122|) . but the parameter dependences are different. Taking account 
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of uncertainties in these lifetime estimates coming from the strong hidden gauge sector, 
we do not regard these estimates to be robust, but we mention possible cosmological 
effects of hidden glueballs if the lifetime happens to be close to 1 sec. In any case, it is 
reasonable to suppose that the dynamical scale of the hidden gauge group should satisfy 



in order for the glueball lifetime to be shorter than 1 sec. Note that there exist many 
excited states of hidden glueballs but they decay into a ground states of scalar or pseudo 
scalar glueballs with lifetime comparable or shorter than the above estimate [39]. 

3.7 Hidden gaugino 

The mass of hidden gauginos are dominantly generated by two-loop RGE running effects 
(see Appendix A), and they are 15-25 GeV at the SUSY scale (~1.5 TeV). However, the 
hidden gaugino mass becomes large when the renormalization scale becomes lower. At 
the one-loop level the following relation holds. 



If we demand Ah ~ 5 GeV, the running mass of the hidden gaugino becomes about 20-50 
GeV. Typical values of the running hidden gaugino masses are shown in Fig. [51 In the 
case where Mh < Mz/2 is satisfied, the Z boson can decay into hidden gauginos at the 
one-loop level. However, the branching ratio is sufficiently small, and the constraint from 
the invisible Z decay width ^30] is easily avoided. 

3.8 Extra Matter 

The typical mass of the extra lepton, m^^, , is 200-400 GeV with k^i = 0.002 (see. TableJ2]). 
Since a large value of ki' is not allowed by the constraints on the vacuum structure (see 
I3.3.ip . the mass of the extra leptons cannot be as large as ITeV. On the other hand, 
the mass of the extra-quark, m^^,, is predicted to be 3-4.5 TeV. The size of ko', and 
hence ttl^^, is limited by the requirement for the successful EWSB and non-existence of 
unbounded vacuum. 



Ah > C(l) GeV, 



(25) 



Mnim^^,) ^ Mh{Mh) 



(26) 
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Q (GeV) 

Figure 5: The running masses of hidden gauginos are shown. The thin sohd hne corre- 
sponds to SU(3)ij case, while the dotted hne corresponds to the case of SU(2)ij. The thick 
sohd hne represents Q = M^, where Q is a renormahzation scale. The gauge coupling of 
S\J{Nh) is taken to be gniQ = l.STeV) = 1.392 for Nh = 2 and 1.122 for Nh = 3, which 
corresponds Ah ~ 5 GeV. 
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As for the scalar states, there are additional contributions from soft SUSY breaking 
masses generated by gauge mediation effects. Since the SM charges of D' are the same as 
those of down- type (s) quark, the soft breaking mass (m^,) is predicted to be ~ 1.5 TeV. 
Similarly, the soft mass of the extra slepton, m^,, is almost the same as that of SU{2)l 
sleptons. 

In most of the viable parameter regions, the lighter mass eigenstate of the extra squark, 
D[ is a bit lighter than its fermionic partner, D', due to the mixing induced by Ao'- The 
extra sleptons are heavier than the extra-leptons, and their masses are about 500 GeV. 

4 Cosmological issues 

In the present model there are lots of cosmological issues to be discussed. Because of the 
hidden gauge symmetry, there appear stable or long-lived particles which may potentially 
spoil success of the standard cosmology. We address these issues in this section. 

4.1 Singlino 

The lightest SUSY particle in the NMSSM sector, except for the gravitino, is the lightest 
neutralino, which mostly consists of the singlino (5) whose mass is around 60 GeV for 
typical interesting parameters. On the other hand , the hidden gaugino {(Jh) also has 
a mass of around 40 GeV for typical parameters (see Sec. 13. 7p . Cosmological scenarios 
significantly change depending on which is the lighter. We discuss cosmology of each case 
separately. 

First, let us assume that the singlino is lighter than the hidden gaugino. The main 
channel of the singlino annihilation is into singlet higgs bosons, ai and hi, both of which 
mostly consist of singlet. However, the annihilation cross section is significantly sup- 
pressed due to the smallness of the coupling k(~ 10~^) in our model (cf. [41j). Therefore 
the singlino decouples from thermal bath soon after the temperature becomes lower than 
the mass of the hidden gauginos. Thus the singlino abundance is huge and it may even 
dominate the Universe soon after the freezeout. The singlino decays into the gravitino 
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and ai with a lifetime of 

If the singhno dominantly decays into the gravitino, a large amount of the gravitino is 
produced. Apparently such large energy density conflicts with the WMAP observation, 
even if there exists a dilution process due to the late-time entropy production. Thus it is 
difficult to reproduce the standard BBN in this setup. 

Next, let us consider the opposite case : the hidden gaugino is lighter than the singlino. 
This is actually the case in the model points PI, P2 and P3 in Tabled In this case, the 
above csomological difficulty can be avoided because the singlino can decay into hidden 
gluino and hidden gluon through the dimension 5 operator!^ 

C= f d^e^^^ + h.c. (28) 

where an = gj^/Air, N = 5 is the number of vector-like matter fields, and W^""^ is the 
gauge supermultiplet of the hidden group. The lifetime is given by 



(29) 

The singlino lifetime is much shorter than the cosmological time scale of the hidden gaug- 
ino decoupling. Thus only the hidden gaugino is left after the freezeout whose abundance 
is sufficiently small and decays into gravitino before BBN begins as described Sec. 14.21 
There is no cosmological difficulty associated with the singlino in this case. Hereafter we 
only consider this case. 

4.2 Hidden glueball and gaugino 

Next, we discuss cosmological consequences of the hidden gauge sector. After the phase 
transition occurring at T ~ Ah, hidden gluons form color singlet states, so-called hidden 
glueballs. As shown in Eqs. (122]) f l24p . the lifetime of the hidden glueball can be short 
enough to be free from the BBN constraint but can also be long enough to modify thermal 
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We neglect the SUSY breaking effects, which does not change the resuh significantly. 
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history before BBN. Hereafter we take the typical scale of the hidden QCD to be Ah ~ 
10 GeV. If it is much smaller, the glueball lifetime becomes too long. If it is much larger, 
the induced CP angle causes a too large EDM (see Sec. 13. 2p . 

Glueballs are formed at the temperature T ~ Kh- Hidden gluons keep thermal equi- 
librium with the SM sector until D' and L' are decoupled at T ~ m^^,/20. This is not 
far from A/^, and hence we assume hidden gluons were in thermal equilibrium at the 
formation of glueballs, although this assumption does not affect the following result at 
all. At the formation, the energy density of thermal plasma is around ~ (7r^5f=K/30)A|^ 
and that of the glueball is around ~ m^^A^ with m^^ being the glueball mass. Since 
mgjj > Khi the hidden glueballs soon begin to dominate the Universe. Its lifetime was es- 
timated in Sec. 13. 61 There we have shown that the lifetime can be shorter than 0.1 sec for 
Aj/ > 0{1) GeV. If the lifetime of the glueball is shorter than 10~* sec, which corresponds 
to a temperature above 10 GeV, it decays as soon as it is formed, and has no effects on 
cosmology. On the other hand, if the lifetime is much longer than 10~® sec but shorter 
than 0.1 sec, the hidden glueball dominates the Universe until it decays at T = T^. In 
this case, the number density of all relics existing before the phase transition, such as the 
gravitino and extra matter to be discussed as well as the baryon asymmetry, are diluted 
by the factor ~ {Tj_/Ah)- We will comment on cosmological effects of late time hidden 
glueball decay in Sec. 14.41 

The hidden gaugino (gn) effectively annihilates into hidden gluon pair at the freezeout 
and the abundance is significantly reduced. After the phase transition, a hidden gaugino 
may form a bound state with a hidden gluon and such a bound state may effectively find 
a partner due to its large geometrical cross section. Then the abundance may be further 
reduced. In any case, the hidden gaugino eventually decays into the gravitino 4'3/2 and 
the hidden gluon. The lifetime is estimated as 

.fe..„.,.3.).I.x.0-^.c(5£|X)^(^)^ (30, 

The gravitino abundance produced by the decay of hidden gluino is sufficiently small. 
Therefore, it does not have significant cosmological effects. 
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4.3 Extra matter 



The lightest particle in the sector of extra matter is stable due to the conservation of the 
hidden gauge charge. The colored particle, D' , is heavier than the lepton-like one, L', 
since koi > k^/. Scalar components obtain SUSY breaking effects from gauge mediation. 
The typical mass relation here is > rri^^, > > ^i'^^ > "^Vl" where rrijj,^ and 
("^L' ^L'^ mass eigenvalues of scalar components of D' and D' ( L' and V). 
We discuss fate of these particles. 

Here we make a brief comment on the splitting of the electrically neutral and charged 
component of L' . The electrically charged component of L', denoted by L^', is slightly 
heavier than its neutral component, L°', due to electroweak radiative corrections, by a 
typical amount of ~ lOOMeV [l2lll3]. Thus L~' can decay into L^' and ez/ with a lifetime 
much shorter than the cosmological timescale. Moreover, its scalar partner, L°', decays 
into LP' by gaugino exchange if L"' is not the lightest i?-odd particle. Since this decay 
process occurs very quickly, the electrically charged component of L' is of no importance 
in the following discussion. Similarly, for the D' multiplet, the lighter scalar component 
D'^ is left among the supermutiplet. 

Now let us follow the evolution of extra matter, L"' and D' . (For notational simplicity, 
we write as L' and D'l as D' hereafter.) It effectively annihilates into hidden gluon 
pairs at the freeze out. By using the cross section of ~ T^ajj / rn"^^^ , the resulting abundance 
is estimated as 



A similar expression holds for D'. Soon after the freezeout, the hidden QCD phase 
transition takes place at T ~ A^^, and then these relic L' and D' particles are connected 
by strings of hidden gauge fluxes. Since there are no other light matter which can form 
pair with L' and D', this matter-string system is stable against creations of L'-L' or D'-D' 
pairs. Thus we need to examine whether this kind of system eventually annihilates or 
not. 

First, we study the dynamics of V-L' pair connected by strings. The mean separation 
i among V particles at the epoch of formation is estimated to be £ ~ ^l/^^ ~ Y^/^'^T"^ ~ 
lO^A^^ The energy stored in one string is thus E ~ A|^4tr ~ 10^ A//, which is consider- 




(31) 
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ably larger than m^^,, with £str ~ ^ being the initial typical string length]^ The energy 
density of the hidden string network is estimated as pstr ~ ^'h 4tr^-' ~ 10-iOA|„ which 
is much smaller than the total energy density of the Universe. String forces attract the 
L'-L' pair, and it easily accelerates the L' to the energy of m^^,). In order for it to 
be annihilated, it must lose its energy so that it can form a bound state with L'. There 
exist several energy loss processes of the matter-string systems during the oscillation : (i) 
hidden gluon emission and (ii) scattering with background plasma. Let us see them in 
detail, 
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(i) Hidden gluon emission 

One of the possible energy loss processes is emissions of hidden gluons. Since the initial 
energy E, carried by L' and L', is much larger than the hidden glueball mass, hidden gluon 
emissions may take place. Note, however, that the acceleration by the string tension is 
around ~ A'jj/m^^,, which is much smaller than the hidden glueball mass ~ Ah- Thus 
emissions occur only when L' and L' meet together on separation of ~ A^;^^. As a rough 
but reasonable estimate, let us assume one hidden glueball is emitted when a L'-L' pair 
passes through each other during one oscillation [38lll5]0 The oscillation period is given 
by, 

tosc ~ I 'J'^^^'^^^H for E = Ajji< m^^, ^g^) 
\ £ for E = Ajji > m^^, 

and here we conservatively assume the emitted hidden glueball energy is ~ Ah- The 

energy loss rate is calculated as dE/dt ~ A/^/tosc- The typical energy loss time is compared 



Initial correlation length of the string (^) is expected to be of order of ^ ~ ^h^- This is much shorter 
than the mean separation of L' particles. Thus, the string connecting them have Brownian motion like 
structure with each step size of ^. Then the string may have length of order of ^str ~ ^^/C 10^*^ Aj^^ 
right after the formation. It might be expected that these small scale irregularities are soon smoothed out 
due to the string tension and string self-interactions, and eventually a straight string connecting L'-L' 
pair is left. Even if we simply regard £str as a maximum and typical initial string length, following energy 
loss mechanisms will dissipate the string energy effectively and the discussion is not affected. 

Gravitational radiations from moving string system also act as an energy loss mechanism, but the 
timescale of gravitational energy loss is much longer than the following processes. See e.g., Ref. [H]. 

Here we assume that the string connects the L'-L' pair in an almost straight line. Otherwise, the 
probability that V finds L' in one oscillation is significantly suppressed. 
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to the Hubble time as 



E/{dE/dt) ^ T^E^"m\ll ^ 10-^3 ^^^V ( ^ V^' ( V'' /lOGeVV 
iHub Mph?jj VlOGeVy VlTeVy VlTeV/ \ Kh ) ' 

(33) 

for E < m^^, , and 

E/jdE/dt) f T V( E y/ lOGeV y 

for E > m^^,. Therefore, the hidden gluon emission effectively dissipates energies of 
strings, if the strings soon relax to a rod-like structure connecting L'-L' pair after the 
formation. 

(ii) Scattering with background particles 

We here note that it is not clear whether a string-matter system relaxes to a rod-hke 
structure so that V and V pass through each other at each oscillation, due to a lack of 
understandings of the dynamics such strings in a cosmological setup. Even if the hidden 
gluon emission does not work, however, there are another process to reduce the energy of 
the system. Since L' {L') has a weak interaction, it scatters off the background plasma 
consisting of leptons and light quarks through tree level interactions. Let us consider 
the interaction of V with background plasma through the Z-boson exchange, when L' 
obtains a large energy E due to the string tension. The typical scattering cross section 



m 



is given by {av) ~ a /{ET) for ET > iri^^, and it loses its energy of order of E 
one collision. Therefore, the energy loss rate is calculated as dE/dt — nT{crv)AE where 
riT ~ denotes the number density of background plasma. By comparing the energy 
loss time scale with the Hubble time, we obtain 

E /{dE/dt) E ,3/ E \ 



Thus the energy of matter-string system is quickly reduced. After the energy decreases 
to ET < rn^^i-i the typical scattering cross section is given by {av) ~ G'^p^ET /m^^,Y and 
an average energy loss per scattering is given by AE ~ TE^/m^ By comparing the 
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energy loss time scale with the Hubble time, we obtain 

Thus the L' would soon be non-relativistic well within one Hubble time. After it becomes 
non-relativistic, the scattering cross section is given by {av) ~ G\T'^ and an average 
energy loss per scattering is I\E ~ TE/m^^, for E < m^^,. Here E is regarded as the 
kinetic energy of L'. Then we obtain 



tHub MpG'pT^ V T y VlTeV. 

Therefore, we conclude that scatterings of L' particles with thermal plasma efficiently 
dissipate the string energy well within one Hubble time. The resultant energy of the 
matter-string system, or the kinetic energy of L'{L'), is expected to be of order of T. 
Since this kinetic energy is comparable or even smaller than aj^m^^,, we reasonably 
expect the formation of L'-L' bound states. Once the bound state forms, the annihilation 
takes place efficiently as shown in Refs. [38l B6] . After all, all the extra matter disappear 
after the hidden QCD phase transition. 

Almost the same arguments hold for D'-D' pair connected by strings. Since D' has 
color and electric charges, its energy loss is expected to be much more efficient than the 
L' studied above. Thus D'-D' pair also annihilates soon after the hidden QCD phase 
transition. 

A complexity arises for D'-L' and D'-L' pairs connected by strings. After they lose 
energies by the same processes described above, they form bound states. They are stable 
because of the conservation of the hidden baryon (lepton) numbers. Since it has a QCD 
color charge, it further forms a bound state with d{d) after the QCD phase transition. The 
bound states D'L'd and D'L'd find each other with the cross section of Aq^j-,. After the 
collision, the system loses energy by the photon emission and finally annihilates rapidly 
[17] . Eventually, the relic abundance of D'L'd and D'L'd bound states is given by [17] 

This is much smaller than the dark matter abundance, but it is still subject to a constraint 
from BBN [IHl HQ]. Taking account of the dilution by the entropy production, this may be 
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marginally allowed. It is also possible to avoid the constraint by introducing the following 
operatoiL^ 

W = ^D'dL'H^, (39) 

with M being the cutoff scale, which may be expected to be O{10^^) GeV. This induces the 
decay of D' into L' with lifetime of 0{1) sec. After the decay, L'-L' bound states connected 
by strings are left with L' particles initially having large kinetic energies. Although the 
energy loss through scatterings with thermal plasma is inefficient at this epoch, they soon 
lose energies through hidden gluon emissions (See Eqs. (j33l) and (l34l)). since the strings 
are straight and they pass through each other in each oscillation. Therefore, after the D' 
decay, the bound states soon disappear. 

4.4 Cosmological scenario 

Now let us summarize the cosmological scenarios in our model. As we have discussed, the 
lifetime of hidden glueball is very sensitive to the model parameters. In the short lifetime 
limit, hidden glueballs do not have significant cosmological effects and we recover standard 
thermal history after reheating. In the long lifetime limit, glueballs once dominate the 
Universe and they finally decay slightly before BBN begins. In this case a significant 
amount of entropy is released. We discuss cosmological scenarios with and without entropy 
production, although the latter is likely to be realized. Let us repeat here that domain 
walls that are formed after electroweak phase transition decay rapidly and harmless in 
our model due to the bias induced by the hidden QCD instanton effect. The singlino 
does not significantly affect cosmology if it is heavier than the hidden gaugino. Extra 
matter particles also quickly disappear after the hidden QCD phase transition due to the 
dynamics of strings connecting them. 

4.4.1 Without entropy production 

First of all, we consider the scenario without late-time entropy production; that is, the 
hidden glueball lifetime is rather short. The reheating temperature, Tr, is bounded above 



This operator can be consistent with the symmetry as weU as the R-parity, with a shghtly 
modified charge assignment: RpiD') = Rp{D') = +, Z^iD') = Zi{L') = 0, and Z:i{D') = Zs{L') = 2. 
Note that this change does not afi^ect the discussion so far at alL 
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from the gravitino overproduction as Tr < 10^ GeV for m3/2 = 0.1 MeV [SO llSTl l52l [531 151] . 
The gravitino can be the main component of dark matter if the reheating temperature is 
close to this upper bound. 

In the presence the PQ sector, the axino (a), a fermionic superpartner of the axion, is 
also produced thermally. The abundance is estimated as [551 ESI ETJ [58] 

This expression is valid also for Tr > 1 TeV. The axino mass sensitively depends on the 
axion model, but it can be as heavy as the gravitino [59j. The axino overproduction is 
avoided for fa > 10^^ GeV if the axino mass is comparable to the gravitino masso The 
saxion (cr), the scalar partner of the axion, is also produced through thermal scatterings 
and coherent oscillation. It generally obtains a mass of at least of the gravitino, and decays 
into two axions with lifetime r ~ 1 x 10^^sec(0.1MeV/mo.)^(/a/10^-'^GeV)^. This decay 
process is constrained from the observation of effective number of neutrino species in the 
cosmic microwave background anisotropy. The reheating temperature of Tr ~ 10^ GeV 
is marginally consistent with this constraint (see e.g., Ref. [60]). But the saxion may 
obtain much larger mass than the gravitino in gauge- mediation [611 [62] . Then the lifetime 
becomes much shorter and the constraint is significantly relaxed. The cosmology of saxion 
depends much on the model of PQ sector, particularly, on how the PQ scalar is stabilized 
at the scale of fa- Here it is sufficient to note that cosmological constraints on the saxion 
can rather easily be avoided. Finally, for the PQ scale of fa ~ 10^^ GeV, the axion 
coherent oscillation is another good dark matter candidate. Thus in this case the dark 
matter may be a mixture of thermally produced gravitino and the axion. 

4.4.2 With entropy production 

A more interesting cosmological scenario appears in the case of late-time entropy pro- 
duction from the hidden glueball decay. Here we assume the lifetime of hidden glueball 
is around 0.1 sec so that their final decay temperature of them is around a few MeV. In 
this case the entropy production dilutes all the abundance of all relics existing before the 



If the axino is lighter than the gravitino, the gravitino can decay into axino and axion, but its hfetime 
is longer than the present age of the Universe. Thus the bound on the Tr does not change. 
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hidden QCD phase transition by up to an amount of A ~ 10'^. Under this circumstance, 
the reheating temperature after inflation can be raised up to, say, Tr ~ lO^^GeV [63] . 
This is because hght gravitinos are thermahzed for such a high reheating temperature 
and their abundance is determined by thermal one, 1^3/2 = 0.417/(7*3/2 where (7*3/2 de- 
notes the relativistic degrees of freedom at the decoupling of the gravitinoo Translating 
into the density parameter, we have VL^/2h^ = 0. 3(7713/2/0. lMeV)(300/(7*3/2)(10^/A). An 
attractive feature is that in such a case, (non-)thermal leptogenesis scenarios [HU |65l 166] 
become viable even under the entropy production |63j . 

Note, however, that if the reheating temperature exceeds the messenger scale, the 
messenger fermions are also efficiently produced and take part in thermal bath. This is 
problematic because the lightest messengers are stable due to a messenger parity !^ One 
can avoid this by demanding Tr < 10^ GeV, where nonthermal leptogenesis still works 
under the entropy production. It is also possible to break the messenger parity to a 
small amount so that messengers can decay. For example, one is allowed to introduce the 
following termq^ 

W = gD,iSd,-^D + gL,iS^^^l, (41) 

with tiny coupling constants go^i and gi^i, and di and £j are the MSSM fields. Even small 
couplings go,! and gi^i of O(10~^°) are sufficient in order for the messengers to decay 
rapidly. 

The axino, as well as the saxion, are also thermalized at the temperature above 
10^ GeV(/a/10^^GeV)^ [HB]- The resulting abundance is similar to that of the gravitino. 
The saxion constraint from the effective number of neutrino species is also marginal, but 
it can easily be avoided by making the PQ scale slightly smaller, or by making the saxion 
heavy as already described. 



Note that ~ 370 in our model if all the matter and gauge fields are thermalized (except for the 
messenger and SUSY breaking sector). 

Also the messenger coupling to the goldstino can enhance the gravitino production [67) . 

Here the Z3 charges of ^ i and q ^ are slightly modified as Zy,{'^ jj i) = 1 and Zz{^ q ,^) = 2, 
respectively. 
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5 Conclusion 



In this paper we have constructed a gauge-mediation model that naturally overcomes 
the ///^//-problem without introducing a large CP angle. The model is based on the Z3 
invariant singlet extension of the MSSM, the so-called NMSSM. The only extension is the 
addition of the extra matter having a charge under a hidden gauge symmetry, as well as 
that of the SM gauge group. In this setup, a desired value of the singlet VEV is induced 
to give a ^-term and the well-known domain wall problem in the NMSSM is solved due 
to the presence of quantum anomaly through the hidden gauge symmetry. This is also 
consistent with the PQ solution to the strong CP problem. Moreover, the relative sign of 
the //-term and the gaugino mass can be preferable for explaining the muon anomalous 
magnetic moment. 

We have also investigated cosmology of the present model. The hidden gauge symme- 
try promises the existence of newly long-lived particles. In particular, the hidden glueball 
may have interesting cosmological implications since it may have a lifetime long enough 
to be a source of late-time entropy production. For the glueballs to decay before BBN, 
the dynamical scale of hidden gauge symmetry cannot be much smaller than 1 GeV. The 
lightest states in the extra matter are also stable, but they efficiently annihilate at the 
hidden QCD phase transition, and are not of cosmological importance. In the case of 
late-time entropy production by hidden glueballs, the reheating temperature can be as 
high as 10^ GeV and (non-)thermal leptogenesis scenario works even after the dilution by 
the entropy production is taken into account, while the gravitino with mass of 0(0. 1) MeV 
accounts for the present dark matter. 

The present scenario has various interesting properties at low energy; weakly coupled 
light scalars and a light neutralino originating from the singlet field, hidden gauge in- 
teraction with a dynamical scale around 10 GeV, and extra vector-like matter at TeV 
scale charged under both standard model and hidden gauge symmetries. These distinct 
features may be tested by the LHC or other experiments in the near future. 
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Figure 6: The corresponding stop mass is Mgusy = ^fn%Q^rrij~ ~ 10 TeV, which shghtly 
depends on the messenger scale. We take A(Msusy) = 0.005, Nh = 3, /c^/ = 0.002, and 
gniMsnsy) = 1-2. The messenger number is taken as = 1. 

Note Added 

If we give up the explanation of the muon g-2 deviation, there is a trivial solution which 
can explain the Higgs boson mass of ~ 125 GeV with ~ 10 TeV stops. In Fig. [6], the 
allowed region of the parameter space is shown for fics < 0, and Agfr = 1000 TeV. Only 
the green region is excluded due to the constraint from the vacuum stability. In the other 
region, the Higgs mass of ~ 125 GeV is explained with ~ 10 TeV stops. There exists an 
almost singlet-like (CP-even) scalar around 300 — 350 GeV, which may be observed at the 
LHC. The solution with /ieg > also exists. But the allowed region is rather small. 
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Appendix A 

Two-loop beta functions for gauge couplings and gauginos are shown. 



dg 



(2)^2 



x=t,b,T 



2 dMj 



dt 



167r2 



x=t,b,T 



where t — InQ, — gn, and 



/ 33/5 + Nh\ 
1 + Nh 
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\ -3Nh + 5 ) 
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9/5 + 3iVH/5 25 + INh 
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6 6 2 

4 4 

\ / 



(44) 



Note that in the case of Nh — 3, 63^^ vanishes accidentally, therefore the terms proportional 
to 63^ and czx are important. We omit contributions proportional to Yukawa couplings 
other than Yt, Yj, and Y^, since they are small. 
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